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The form factors and the coupling constant of the B^BK vertex are calculated using the QCD sum 
rules method. Three point correlation functions are computed considering both K and B mesons 
off-shell and, after an extrapolation of the QCDSR results, we obtain the coupling constant of the 
vertex. We study the uncertainties in our result by calculating a third form factor obtained when the 
B* s is the off-shell meson, considering other acceptable structures and computing the variations of the 
sum rules' parameters. The form factors obtained have different behaviors but their simultaneous 
extrapolations reach to the same value of the coupling constant Qb%bk = 10.6 ± 1.7. We compare 
our result with other theoretical estimates. 
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I. INTRODUCTION 

Hadronic three-meson vertex functions or form factors are basic inputs to phenomenological theories of nuclear 
processes. They play an important role in cross section calculation, which may vary in a few orders of magnitude 
if they are considered. The form factor depends on the momentum Q 2 and some parameters. As usual, the form 
factor is parametrized by ad hoc function, which is not unique and can be Gaussian, exponential or monopolar. From 
the parametrization, we can extract a cutoff parameter A which is associated with the sharpness of the form factor. 
Therefore, knowing its precise functional form and its cutoff parameter are essential. 

Some years ago, we started a program to compute form factors with the QCD sum rules (QCDSR) Q. Our group 
has developed a method which leads to less ambiguities in the determination of the coupling constant. The method 
consists of computing two form factors of the same three meson process: one form factor is evaluated when the heavy 
meson is off-shell and the other form factor is obtained when the light meson is off-shell. By extrapolating both form 
factors to the momentum Q 2 equal to — m^, where m 2 ^ is the mass of the off-shell particle, we obtain the coupling 
constant of the vertex. After our pioneer DDp work 0, we concluded that when the form factor tends to be harder as 
a function of Q 2 ,the off-shell meson in the vertex is heavy, which means the cut-off parameter is larger. A persistent 
study of vertices involving charmed mesons (D*Dtt H£|, DDJ/ip [3, D*DJ/ij) @], D*D*tt @, Q, D*D*J/ij) (§], 
D S D*K, D*DK and D*D*p [H|), motivated by the development of effective theories related to charm particles 
interaction, showed a similar conclusion in all of them. 

However, in recent years, due to the precise measurements of B decays performed by BELLE, BES and BABAR, the 
physics of B meson has gained a new relevance. The observed -B s i(5830) by CDF Collaboration [l2| and S S 2(5840) by 
CDF and DO Collaboration [l3| stimulated our interest in the physics of bottom mesons interactions [l4j]. In particular, 
the unobserved B* Q (5725) state and the axial B Sl (5778) have been interpreted as hadronic bound states (hadronic 
molecules), because their masses are close to the thresholds of the corresponding hadronic pairs. In reference [l5j |. 
the authors used a phenomenological Lagrangian approach to calculate the strong and radiative decay of these new 
mesons, _B* o (5725) and £? Sl (5778), which were considered as bound states of the BK and B*K mesons respectively, 
to provide more information about some properties. In particular, the coupling constant of B*BK vertex was an 
input in this calculation. 

Considering that the form factors are introduced in the development of effective theories to study bottom interac- 
tions, in this paper, we calculated the form factors and coupling constant of the B*BK vertex, using the QCDSR 
method, which is the only method that permits to extract the coupling constant of the three meson processes without 
any dependence on other empirical coupling constant. Therefore, we evaluate two form factors, one when B is the 
off-shell particle and another when K is the off-shell and extrapolate these results to obtain the coupling constant 
9b*bk- We also want to estimate the uncertainties of the QCDSR method by analyzing different sources of errors. 

In section II, we describe the QCDSR technique; in II. A, the QCD side for this vertex and in II. B, the phenomeno- 
logical side. In section III, we show the results; in section IV, we estimate the uncertainties and, finally, conclude. 



II. THE SUM RULE FOR THE B*BK VERTEX 

Following our previous works, specifically in Ref.[l(j, the three-point function associated with the B*BK vertex is 
givem by 

T^W) = / d^xd^y e^-*e-^'^-y(0\T{j?(x)j^(y)jP(0)}\0) (1) 

for an off-shell B meson, and: 

r<*W) = J d^xdS e**> e-^-v)-v ( \T{j^(x)jP(y)j^(0)}\0) , (2) 

for an off-shell K meson. The expressions for the vertices ([I} and ([2]) contain different numbers of Lorentz structures, 
and each structure can be a different Sum Rule (SR) to be considered. Equations ([T]) and ([2|) can be calculated in 
two different ways: using quark degrees of freedom -the theoretical or the QCD side- and using hadronic degrees 
of freedom -the phenomenological side. In the QCD side, the correlators are evaluated using the Wilson operator 
product expansion (OPE), which incorporates the effects of the QCD vacuum through an infinite series of condensates 
of increasing dimension. The phenomenological side is written in terms of hadronic degrees of freedom, which is the 
responsible for the introduction of the form factors, decay constants and masses. After performing a double Borcl 
transformation, both representations are matched invoking the quark-hadron global duality. 
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The OPE side 



In the OPE or theoretical side, each meson interpolating field appearing in Eqs. ([TJ and (|2|) can be written in terms 
of the quark field operators. For the B off-shell case, Eq. ([T]), we use the following meson currents: 

3 (°) = h^s 



J 



(3) 



and when the K is off-shell, we have chosen the pseudoscalar current for it, following the same procedure as in previous 
works as D*Dtt [| and D* S DK [T3]: 

jK{x)=iq lb s, (4) 

this choice of current is motivated by the coupling constant that is obtained by extrapolating the <?s*sif (*? 2 ) ^° 
momentum Q 2 equal to the square mass of the particle, that in this case will be rn 2 K . Otherwise, if we use the axial 
vector current when K meson is off-shell, the sum rule obtained is also good and the coupling constant result is very 
similar to the pseudoscalar case. 

The u, d, s and b are the up, down, strange and bottom quark fields, respectively. Each one of these currents has 
the same quantum numbers of the associated mesons. 

In order to obtain the QCD side, we use the Cutkosky rule and we obtain an expression where all the possible 
structures appear. We can obtain a sum rule for each structure that have an invariant amplitude, which can be 
written in a double dispersion relation over the virtualities p 2 and p' 2 , holding Q 2 = 



fixed: 



-i poo />oo 



(III 



p T {s,u,Q 2 ) 
(s-p 2 )(u-p> 2 ) ' 



(5) 



where p T (s,u,Q 2 ) equals the double discontinuity of the amplitude T T (p 2 , p' 2 , Q 2 ) , T index means B off-shell or K 



off-shell, q = p' — p and s m in and \ 
for B off-shell and s m in = m 2 and 



are integration limits equal to s r 
= t 



ml for K off-shell. 



(nib + rn s ) and 



t 



The invariant amplitudes receive contributions from all terms in the OPE. The first of those contributions comes 
from the pcrturbativc term, which is represented in Fig. [T] 



Bf 



<ln 



Pv 



K~ 



B+ 



K~ 




Pv 



Bf 



FIG. 1. Perturbative diagrams for the B off-shell (left) and K off-shell (right) correlators. 



In order to obtain the form factor, we choose one of the structures in Eqs.(Q} and ([2]), the one that has less 
ambiguities in the QCD sum rule approach, which means less influence from the higher dimension condensates, better 
stability as a function of the Borel mass and has larger pole than continuum contribution. The other structures that 
are "good " sum rules can also be considered to improve the uncertainties of the method. For B off-shell meson, 
which has two different structures, we choose the p'^p'^ because it is the best structure with less ambiguities. For K 
off-shell meson, we have p v and p' v structures, both being excellent sum rules. The chosen is p' v and the other will be 
used in the calculation of the uncertainties. 

The corresponding perturbative spectral densities, which enter in Eq. ([S]), arc: 



p( B \s,u,Q 2 ) = 
for the p'uP'v structure of the B off-shell case, and 



[(2m b (A-B))) 



\s,u,Q 2 ) = - 



B 



m b m s 



2m 2 



(0) 



(7) 
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for the p' v structure of the K off-shell case, where A = A(s, u, t) = s 2 + t 2 + u 2 — 2st — 2su — 2tu, s = p 2 , u = p' 2 , 
t = —Q 2 and A and B are functions of (s, u, t), givem by the following expressions: 

A = -^L(3cos0-l); S = 27rHcos0; (8) 



where 



J\ = ml cose = _Hfco-"_j~W. 

m 1*1 

with i = s, Tf = and e = 1 for B off-shell and i = b, r\ = — 1 and e = for K off-shell. 

The main OPE contribution is the perturbative. The next terms in the OPE expansion are the quark condensate, 
the gluon condensate and the mixed quark-gluon condensate. In this work, the quark condensate was calculated and 
has the following expression obtained for the K off-shell case: 

T< bb> = m s <bb> e m * /M V (10) 

this contribution is proportional to the bottom condensate, consequently, the contribution is null The other two 
quark condensates contributions for the K off-shell case vanish after the application of the double Borel transformation. 
For B off-shell case, the quark condensate does not contribute for the structure considered. In this calculation, we 
do not include gluon condensate and mixed quark gluon condensate, because as seen in previous calculation of the 
form factors of D*Dtt and B*Bn vertices [||, we found out that the gluon contribution is negligible as compared with 
the perturbative one. Moreover, the gluon condensate decreases with the Borel mass and the most important role 
of the gluon condensate is to guarantee the stability as a function of the Borel mass. In this vertex, this stability is 
guaranteed together with a good contribution of the pole. 



B. The phenomenological side 

The three-point functions from Eqs. {TJ and ([2]), when written in terms of hadron masses, decay constants and 
form factors, are the phenomenological side of the sum rule, which is based on the interactions at the hadronic level. 
These interactions are described here by the following effective Lagrangian [l6j : 



C 



B*BK = 19B* BK 



B^iBO^K - d^BK) - B^iKd^B - d^BK) 



[11) 



from where we can extract the matrix element associated with the B*BK vertex. 

The meson decay constants fx, fg and fg* are defined by the following matrix elements: 



{0\j?\K(p))=if KPv , 



(12) 



(0|jf \B* s (p)) =m B: fg : e*Jp) 



(13) 



and 



(14) 



where e* is the polarization vector for the B* meson. Saturating Eqs. ([IJ and <j2j) with K and B* states and using 
Eqs. CE2]) - (fig]), we arrive at 



Tl B J(p,p / , q )=g < B B J BK ( q 

Ip'uPv - p'^Pv 1 + 



fBifufBrng* 



{p 2 - m 2 B: ){q 2 - m 2 nB ){p" - m 2 K ) 
(m% — q 2 ) , 

1 continuum 



(15) 
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when B is off-shell. Using the matrix element of K meson equal to 



(0\j?\K(p)) = '-^f K , 



(16) 



we arrive at an expression for K off-shell: 

Ti K \p,p' >q )=-g( B % K (q 2 ] 



fBfKfB*rn 2 B ,m 2 B m 2 K 



m b m s (p 2 - m 2 B )(p' - m 2 B ,)(q 2 - m 2 K ) 
(m| - q 2 ) 



m% 



"continuum". 



(17) 



In the Eqs. (p~5|) and dTTJ) , we can see all the structures that appear in the QCDSR and the "continuum" contribution, 
which may be subtracted assuming that its contribution is equal to the contribution of the QCD side, by introducing 
the continuum thresholds so and uq [l7| as upper integration limits in Eq. ([5]). 

In order to improve the matching between the two sides of the sum rule, we perform a double Borcl transformation 
[ijj in the variables P 2 = —p 2 — » M 2 and P' 2 = -p' 2 -> M' 2 , on both invariant amplitudes r' T ' (QCD side) and 



(phenomenological side). Equating the results, we obtain the form factors g B J B x(Q 2 ) that appear in Eqs. ([15 
and (JT7J). 



(T) 



III. RESULTS AND DISCUSSION 



Table U shows the parameters used in the present calculation. We have used the experimental value for fx of 
Ref . [lU , for /s» and fs from Ref . [l!| ; for m s from Ref . [2(j and for mf, we used recent results of Ref . [HI, [22| . The 
continuum thresholds are givem by sq = (m, + A s ) 2 and uq = (m + A u ) 2 , where m, is the mass of the incoming 
meson, m is the mass of the out coming meson and A u and A s are usually around 0.5 GeV, but this value can be 
varied to obtain a good pole-continuum contribution and stability of the sum rule with the Borel mass parameter. 





q 


s 


b 


K 


B 


b; 


m (GeV) 


0.007 


0.13 


4.20 


0.49 


5.27 


5.41 


/ (MeV) 








160 


208 


250 



TABLE I. Parameters used. 



A. B off-shell form factor 



In order to obtain the B off-shell form factor, we choose the structure p'^p'^ which has an excellent stability and 
its pole contribution is larger than the continuum contribution in a Borel region between 10 < M 2 < 24 GeV 2 . The 
first order in the OPE, the quark condensate does not contribute in this structure. In this paper, the next order in 
the OPE, the gluon condensate was not calculated, because in our previous work of D*Dn form factor [3| and also 
in J/^D*D [23], we found out that the gluon condensate contribution is less than 10% for a Borel mass of 5 GeV 2 . 
In this case, we are dealing with heavy mesons and so we arc working with a Borel mass that is 21 GeV 2 with 
which wc arc expecting the gluon condensate to be negligible and it does not contributing to modify the result of the 
pcrturbative term. 

Performing the match beetwen the QCD and the phenomenological side, in Fig. we show the Borel window 

stability for different value's combinations of thresholds and in the same Fig. (Hb), we can see the pole and continuum 
contributions observing that the pole contribution is always bigger than the continuum contribution in the same Borcl 
window of stability, therefore we have good credibility for this sum rule. Joining this criteria with the thresholds 
A s = 0.6 GeV and A. u = 0.5 GeV, we obtain the form factor for B off-shell: g^BKiQ 2 )' where we have used 

2 2 

■jtfrx = — 2— r — 2" as a Borel mass relation, which is usual when light and heavy mesons are involved. 

The triangles in Fig. © correspond to the sum rule result, which can be fitted by the monopolar function represented 
by the dashed line in Fig. ([5]) and given by 
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■ as=au=0.5 

■ AS=AU=0.6 
-As=0.6; Au=0.5 



a Pole 

• Continuum 







M 2 (GeV z ) 
(a) 



M (GeV ) 
(b) 



FIG. 2. a) g { B* B AQ 2 = 1 GeV 2 ) as a function of the Borel mass M 2 and different thresholds and b) pole-continuum contribu- 
tions. 



The coupling constant is obtained when the form factor is extrapolated to Q 2 = — mf., where tot is the mass of the 
off-shell meson @, [1(1 EH- Using Q 2 = —m 2 B in Eq (fT8|) , the resulting coupling constant is equal to: 

9b : bk = 10.6 (19) 

For the B meson off-shell case, we have other structure, p^p^, to work with. This sum rule has quark condensate 
contribution but it has not good stability and its continuum contribution is larger than the pole contribution, therefore 
this sum rule is not considered a good sum rule. 



B. K off-shell form factor 



In the case K off-shell, we have two structures in Eq. (fTT|) that can be used: p' v and p v . Both structures give good 
sum rule results, that means a good pole-continuum contribution and good stability. The quark condensate of the 
bottom is considered, but the contribution is null. At first, we have chosen the p' v structure to obtain the form factor. 
The stability for different thresholds are showed in Fig. ([3k,). For A s = 0.7 GeV, A M = 0.7 GeV and for a Borel mass 
of 7.2 GeV 2 , we obtain a pole contribution of about 70% of the total contribution ( Fig. ([3]) part b)). 



■4S=AU=0.5 
■ AS^AU-0.6 
AS=AU=0.7 



M (GeV ) 

(a) 




M (GeV ) 
(b) 



FIG. 3. a) g ( B» ] BK {Q 2 = 2 GeV 2 ) stability as a function of the Borel mass for different threshold values, structure p v and b) 
pole and continuum contributions. 
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The white circles in Fig. (j6)) represent the sum rule result fitted by a Gaussian function that is the short dashed 
line in the same figure and is given by 



by extrapolating Eq. ([20)) to Q 2 = —m 2 K , we obtain the coupling constant: 

9b%bk = 10.7 



(20) 



(21) 



We can see an excellent agreement between the results obtained in Eqs. ([191 and (|2fj) . as expected. 

We noted the same vertex behavior concerning the shape of the form factor that we obtained in previous works, 
when the heavy meson is the off-shell particle, the form factor is "harder" and has a big cutoff parameter A, which is 
defined from Eqs. (| 18(1 and (|20| . which have a general Gaussian form: 



9B'Bi< 



(g 2 ) = A e -w 2 /A) 2 



(22) 



and a monopolar form 



9b*bk{Q 2 ) 



A + Q 2 

where A parameter is equal to 34.99 GeV 2 for B off-shell meson and is 7.49 GeV 2 for K off-shell meson. 



(23) 



IV. UNCERTAINTIES OF THE QCDSR METHOD 

The uncertainties in the QCDSR calculation come from different sources such as decay constant errors, mass 
uncertainties, other "good " structures which also can be considered and the QCDSR variation parameters themselves. 
In the next section, we discuss each case separately in sub-sections and use the results to find an error of the method 
in our result. 



A. B* off-shell 

At first, we calculate a third sum rule with the other heavy meson of the vertex off-shell: the B*. Since this meson 
has a small mass difference when compared with the B meson mass, we expect to obtain a similar result to the one 
obtained for B off-shell. We have two structures, p'^p' v and p^p'^, to work with in the method. The p^p'^ structure 
does not present good stability in the Borel mass, therefore it is not considered. The other structure, p'^p'^, has good 
stability and the pole is larger than the continuum contribution, for A s — 0.7 GeV and A u = 0.4 GeV. Fig. |4]a) 
shows the Borel mass stability and in b) the pole-continuum contribution. 



■ AS=AU=0.4 
■AS=AU=0.5 

AS=AU=0.6 

AS=AU=0.7 

-As=0.7; Au=0.4 



M (GeV ) 
(a) 




10 12 14 16 18 20 22 24 



M (GeV ) 
(b) 



FIG. 4. (a) Stability of g^J^iQ 2 = 1 GeV 2 ) as Borel mass and for different thresholds and (b) pole-continuum contributions. 
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We extrapolate the QCDSR results, the dots in Fig. |6l by using a monopolar function until reaching Q 2 = —m 2 B » , 
and wc obtain 



93.03 



37.80 + g 2 



for the form factor, and 



9 b* bk = 10.25 



(24) 



(25) 



for the coupling constant. Both results, Eq. (|24| and Eq. (|18| show the same extrapolation function with similar 
parameters. 



B. K off-shell, structure p v 

For the light meson off-shell, K, we also have the p u structure, which gives a good sum rule and is used here 
to calculate the uncertainties of the method. For this structure, we analyze the Borel mass stability for different 
threshold parameters, Fig. [5] a), and the pole-continuum contribution, showed in Fig. [5]b). Wc obtain a Gaussian 





(a) 



(5 b) 



FIG. 5. a) Stability of g { gJ BK (Q 2 = 2 GeV 2 ) as a function of the Borel mass, for different threshold values and structure p„; 
b) Pole and continuum contributions. 



+f? BK {Q 2 ) = 10.65e-W 2 / 4 ") 2 , 



(26) 



function by extrapolating the QCDSR result given by: 

9b*bk( 

and making Q 2 = —m 2 K , we obtain the coupling constant of the vertex which is equal to: 

9b*bk = 10.7 (27) 
In Fig. ([6]) we summarize the results obtained with both structures showing an excellent agreement. 



C. Uncertainties of parameter variations. 

In order to know, how much docs the QCDSR depend on the parameters' variations, we analyze, in this work, the 
experimental data uncertainties for the decay constant Jk [H| , the theoretical uncertainties for m s , m,b [HI HH , 
}b and Jb* [H| and the variations over the QCDSR parameters that are the Borel mass, M 2 , and the thresholds As 
and Au. 

To compute the error, we proceed in the following way: after calculating the coupling constant, we compute a new 
coupling constant, but now with all parameters kept fixed, except one, which is changed according with its intrinsic 
error, given in Table [TTJ Then, we move to the next parameter to be varied, keeping all others fixed. At the end of the 
procedure, we know how sensitive is this vertex respect to each parameter. In table [ill that contain the percentual 
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Deviation % 


Parameters 


B off-shell 


K off-shell 


B* off-shell 


}k = 159.8 ± 1.4 ± 0.44 (MeV) 


1.04 


1.04 


1.02 


fa = 208 ± 10 ± 29 (MeV) 


15.92 


15.97 


15.94 


f B * = 250 ± 10 ± 35 (MeV) 


15.22 


15.25 


15.29 


M 2 ± 10% (GeV) 


6.05 


5.9 


2.41 


m b = 4.20 + 0.17 - 0.07 (GeV) 


13.69 


28.80 


14.18 


m s = 104 + 26 - 34 (MeV) 


8.98 


14.67 


5.59 


As ±0.1 e Au±0.1(GeV) 


13.51 


2.79 


14.30 



TABLE II. Percentage deviation related with each parameter. 



deviation for the three cases (B, K and B* off-shell mesons), we observe that the biggest variation on the coupling 
constant is produced when fs, /s» and mb were varied. Finally, the uncertainties for each coupling constant, which 
is obtained with B, K and B* off-shell are: 

9b:bk = 10 - 6 ±1-5 (28) 

for B off-shell meson, 

g%) BK = 10.3 ±1.7 (29) 

for K off-shell, where were used both structures that give good SR, and 

9%i K = 11.0 ±1.5 (30) 

for B* off-shell. Fig.© summarize these results, where the error bar represents the uncertainties for each calculation. 

Considering the three off-shell cases, for all possible and "good" sum rules structures, we obtain the coupling 
constant of the B*BK vertex. After computing all the uncertainties, we obtain the mean value: 

g B * a BK = 10.6 ±1.7 (31) 




FIG. 6. The three cases considered here: off-shell K, B, B* and for all possible sum rules. 
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The uncertainties of the QCDSR coming from different sources (other structures with the same mesons of shell, 
other heavy meson off-shell, the extrapolation fit of the sum rule result and the errors in masses, decay constants, 
condensates, the choice of the Borel mass and the continuum threshold parameters) is near 16 %. The usual errors 
in the QCDSR are around 20% [H[. 

Also, from Eqs. (fT8|) , ((20)) and ([24]) , we can extract the cut-off parameters, A, associated with the form factors. 
In table IIII[ the cut-off parameters are showed and we can observe that they follow the same trend as we had been 
observed in Refs.JH, @: the value of the cut-off is directly associated with the mass of the off-shell meson probing the 
vertex. A harder form factor means bigger cut-off parameter. 



Meson off-shell 


A parameter (GeV 2 ) 


B 


34.99 


b: 


37.80 


K 


7.49 



TABLE III. Values of the cut-off parameter A. 



Recently, our extrapolation procedure 0- 25| was used to calculate the same coupling constant, gs*BK, by the 



authors of Rcf. [26| . Even if the procedure to obtain the coupling constant was the same, the form factors obtained were 
different. This calculation do not include the analysis of the pole-continuum contributions, exponential extrapolations 
functions with three parameters were used to fit the form factors and different values for the parameters were used in 
the SR, resulting in a different value of the coupling constant. 

There are some possibilities that could be considerer to compare our result with other theoretical predictions. If 
we use arguments of heavy hadron chiral perturbation theory (HHChPT), the couplings for the bottom-light vertex 
9b*bk are related to the charm-light vertex 9d* dk through the relation [27|-[28|: 

m B 

9 b* bk — 9d*dk , 

m D 

where if it is used tub = 5.279, ran = 1.8693 and 9d*dk = 2.84( which is our QCDSR result 0), we obtain 
9b* bk = 8.02, that is in complete agreement with the result of this paper. 

Also, 9b* bk is related to 9b*b-k by 5/7(3) symmetry and the coupling constants became equal 9b* s bk = 9b*Btt- 
In table Hvl we present other theoretical estimatives from QCDSR, Light Cone SR and other methods of the 9b*b-k- 
We find out that SU(3) symmetry present different degrees of violation where is compare with our estimative of the 
coupling constant. The third column in the table llVl show these results and we obtained similar conclusion in our 
review 12511 . 



Approach 


9B* BK 


violation 


QCDSR[30] 


20 ±4 


50% 


QCDSR[30] 


15 


44% 


LCSR [31] 


28 ±6 


65% 


QCDSR [32] 


14 ±4 


25% 


LCSR [33] 


22 ±9 


55% 


QCDSR [34] 


42.5 ±2.6 


76% 


QCDSR plus meson loops [35] 


44.7 ± 1.0 


73% 


dispersive quark model [361 


32 ±5 


69% 


Dyson-Schwinger equations [371 


30.0±?;4 


73% 



TABLE IV. Summary of theoretical estimates for gB'Bw 



Concluding, we calculate the form factors of the vertex showing that they have two different forms, one when the 
heavy meson is the off-shell particle and another when the light meson is the off-shell one, with their extrapolations 
to the on-shell point leading to the coupling constant. 

In this paper, we made all the calculations taking control of the SR stability and the pole contribution is always 
contributing more than the continuum. Besides, we computed all the "good " sum rules and we obtained the same 
value of the coupling constant. Furthermore, we computed the sum rules taking into account the errors in masses, 
decay constants, condensates, choice of the Borel mass, continuum threshold parameters, obtaining a final result with 
errors that are near 16%. 
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